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ABSTRACT 

Objectives: To characterise the mRNA expression 
patterns of early and advanced stage colorectal 
adenocarcinomas of Malaysian patients. 
Design: Comparative expression analysis. 
Setting and participants: We performed a 
combination of annealing control primer (ACP)-based 
PGR and reverse transcription-quantitative real-time 
PGR for the identification of differentially expressed 
genes (DEGs) associated with early and advanced 
stage primary colorectal tumours. We recruited four 
paired samples from patients with colorectal cancer 
(GRG) of Dukes' A and B for the preliminary differential 
expression study, and a total of 27 paired samples, 
ranging from GRG stages I to IV, for subsequent 
confirmatory test. The tumouric samples were obtained 
from the patients with GRG undergoing curative 
surgical resection without preoperative 
chemoradiotherapy. The recruited patients with GRG 
were newly diagnosed with GRG, and were not 
associated with any hereditary syndromes, previously 
diagnosed cancer or positive family history of GRG. 
The paired non-cancerous tissue specimens were 
excised from macroscopically normal colonic mucosa 
distally located from the colorectal tumours. 
Primary and secondary outcome measures: The 
differential mRNA expression patterns of early and 
advanced stage colorectal adenocarcinomas compared 
with macroscopically normal colonic mucosa were 
characterised by AGP-based PGR and reverse 
transcription-quantitative real-time PGR. 
Results: The RPL35, /?PS25and TIMP1 genes were 
found to be overexpressed in both early and advanced 
stage colorectal adenocarcinomas (p<0.05). However, 
the ARPC2Qene was significantly underexpressed in 
early colorectal adenocarcinomas, while the advanced 
stage primary colorectal tumours exhibited an additional 
overexpression of the C6orf173Qene (p<0.05). 
Conclusions: We characterised two distinctive gene 
expression patterns to aid in the stratification of 
primary colorectal neoplasms among Malaysian 
patients with GRG. Further work can be done to assess 
and compare the mRNA expression levels of these 
identified DEGs between each GRG stage group, 
stages l-IV. 



Strengths and limitations of this study 



■ This regional-based study has a relatively small 
sample size due to the lack of a designated 

k Tissue Bank and colorectal cancer (GRG) patient 
" volunteers. 

■ All subjects were newly diagnosed with GRG, 
and were not associated with any hereditary syn- 
dromes, previously diagnosed cancer or positive 
family history of GRG. 

■ The findings of this study are considered 
preliminary. 



INTRODUCTION 

Cancer staging is vital for patient manage- 
ment, especially in prognosis prediction and 
planning of treatment intervention/ This is 
especially so in the colorectal cancer (CRC) 
staging system. As such, there have been 
many noteworthy improvements since the 
introduction of the classical Dukes' staging 
system, followed by the modified 
Astler-CoUer staging system; to the latest 
seventh edition of tumour-node-metastasis 
(TNM) staging system published by the 
American Joint Committee on Cancer 
(AJCC).^"^ The TNM staging system allows 
the incorporation of various clinical informa- 
tion (which are obtained through histopatho- 
logical examination, radiological imaging 
and surgical findings), for accurate CRC 
stratification.^ However, these clinical assess- 
ments are greatly dependent on the expert- 
ise of pathologists, radiologists and clinicians. 

The TNM classification is applicable for 
both clinical (cTNM) and pathological 
(pTNM) staging of primary colorectal 
tumours. Typically, it involves the assessment 
on the depth of bowel wall invasion at the 
time of diagnosis and the presence of 
regional lymph nodes metastases, as well as 
the presence of distant organ metastasis.^ 
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As a potentially worse patient outcome with more 
advanced disease stage is the core concept in cancer 
staging, AJCC revises the TNM classification system 
every few years with an attempt to formulate it for more 
accurate patient prognostication.^ The latest seventh 
edition has further detailed the subclassification of the 
pN category and the assessment of discontinuous/ satel- 
lite tumour foci. However, these revisions have increased 
the complexity and subjectivity during evaluation, and 
thus might lead to interobserver variability and hamper 
its efficiency in routine clinical practise.^ ^ In addition, 
current clinicopathological parameters are insufficient 
to address the great biological and genetic heterogen- 
eity of CRC in patients' outcome and treatment 
response prediction. From the perspective of clinical 
oncology, the integration of molecular biomarkers into 
existing clinicopathological assessment will further 
refine the cancer management in future. 

Over the past decades, many researchers have 
attempted to establish gene expression signatures spe- 
cifically for the diagnosis, prognostication and recur- 
rence prediction of sporadic CRC. Transcriptional 
profiling promises a fairly dynamic view on the cellular 
functions, regulatory mechanisms and biochemical 
pathways involved in the disease pathogenesis and pro- 
gression.'^ Various gene expression profiling techniques 
ranging from differential display, SAGE to microarrays, 
have been utilised. Despite its wide application in gene 
expression profiling, microarray experiments have 
been subjected to various sources of variability, false- 
positives as well as statistical and bioinformatic chal- 
lenges. To date, none of the molecular markers 
described has been validated and employed in routine 
clinical practise owing to the poor reproducibility of 
the identified differentially expressed genes (DEGs) 
between different profiling platforms.^ Although the 
KRAS mutation and mismatch repair status have shown 
promising prognostic and predictive values, they have 
yet to be incorporated into either routine pathological 
reporting systems or TNM staging systems.^ 

As most of the molecular studies on CRC were based 
in Western populations and different molecular 
changes were thought to underlie the development of 
sporadic CRC in populations with different genetic 
backgrounds, we aimed to investigate the changes in 
mRNA expression patterns in primary sporadic colorec- 
tal tumours with regard to our Malaysian patients. In 
our study, we employed a combined approach of a 
two-step ACP-based PGR and real-time reverse tran- 
scription PGR to characterise the gene expression pat- 
terns for early and advanced stage sporadic colorectal 
adenocarcinomas. 



MATERIALS AND METHODS 

Patient selection and specimen collection 

All patients presented with histologically confirmed colo- 
rectal adenocarcinomas and were staged accordingly to 



Table 1 Cancer staging of recruited subjects 
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the AJCC TNM staging system (table 1). The staging of 
cancer was performed by taking into consideration their 
histopathological reports, CT images, morphological 
evaluations during surgery and serum carcinoembryonic 
antigen levels. All patients were newly diagnosed with 
CRC, and were not associated with any hereditary syn- 
dromes, previously diagnosed cancer or positive family 
history of CRC. Initially, four patients with CRC stages 
I-III were recruited for the preliminary ACP-based PCR 
analysis, while another 27 patients with CRC stages I-IV 
were recruited for subsequent reverse transcription- 
quantitative real-time PCR (RT-qPCR) analysis. The 
patients' group was composed of the three main ethnic 
groups in the Malaysian population, that is, Chinese, 
Malays and Indians, in order to ensure a better repre- 
sentative of the study population. 

The subjects were admitted to the University Malaya 
Medical Centre, Kuala Lumpur, Malaysia, and under- 
went curative surgical resection between 2010 and 2011. 
None had received preoperative chemoradiotherapy. 
The study protocol was approved by the Ethics 
Committee Board of University Malaya Medical Centre 
(Ref. No. 654.1), and written informed consent was 
obtained from all study patients. The tumouric speci- 
mens were excised from the primary colorectal tumours, 
while the non-cancerous tissue specimens were obtained 
from distally located macroscopically normal colonic 
mucosa. Both colorectal tumour and paired non- 
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cancerous tissue specimens were immersed in RNAlater 
RNA Stabilisation Reagent (Qiagen) immediately after 
excision and stored at -80°C. 

Total RNA extraction 

Total RNA was extracted from homogenised colonic 
tissues with the RNeasy Plus Mini Kit (Qiagen) accord- 
ing to manufacturer's instructions. Subsequently, the 
RNA yield and integrity were ascertained via Agilent 
2100 Bioanalyser in conjunction with Agilent RNA 6000 
Nano Kits (Agilent Technologies). The values of RIN 
were then determined in order to assess the integrity of 
the isolated total RNA. In this study, only RNA samples 
with RIN values of 8.0-10.0 and rRNA ratios (28S/18S) 
of 1.5-2.5 were selected for successive applications. 

ACP-based PGR analysis 

A. First-strand cDNA synthesis 

The synthesis of first-strand cDNA was performed 
according to the manufacturer's protocol for the 
GeneFishing DEG Premix Kit (Seegene), as follows: 3 ]ig 
of total RNA was added with 2 viL of 10 \iM dT-ACPl 

(5'-CTGTGAATGCTGCGACTACGATXXXXX (T) 18-30 
and RNase-free water to a final volume of 9.5 \iL. The 
mixture was then incubated at 80° C for 3 min, followed 
by chilling on ice for another 2 min. Subsequently, 4 \ih 
of 5X RT buffer (Mbiotech), 5 >iL of 2 mM dNTP 
(Fermentas), 0.5 ]jiL of 40 U/]iL RNase inhibitor 
(Mbiotech) and 1 \iL of 200 U/\iL M-MLV reverse tran- 
scriptase (Mbiotech) were added. This mixture was then 
incubated at 42°C for 90 min, heated at 94°C for 
another 2 min and chilled on ice for 2 min. Finally, 
80 \il of DNase-free water was added to dilute the synthe- 
sised cDNA. The first-strand cDNA was stored under 
— 20°C until further analysis. 

B. ACP-based GeneFishing PGR 

First, all four cDNA samples within each GRG and 
control group samples were pooled together in equal 
amounts. The characterisation of DEGs was then con- 
ducted via ACP-based PGR based on 20 arbitrary AGP 
primers (Gat. No. K1021) in a thermal cycler 
(Mastercycler Gradient, Eppendorf) according to the 
manufacturer's protocol (GeneFishing DEG Premix Kit, 
Seegene). Initially, the synthesis of second-strand cDNA 
was commenced in a one-cycle first-stage PGR: 94° G for 
5 min, 50°G for 3 min and 72°G for 1 min. Next, the 
constructed second-strand cDNA was subjected to 
second-stage PGR with 40 cycles of a denaturing step at 
94°G for 40 s, annealing step at 65°G for 40 s and exten- 
sion step at 72 °G for 40 s. Finally, a final extension step 
at 72 °G for 5 min was carried out. The amplified pro- 
ducts were then separated on 3% (w/v) agarose gels 
stained with ethidium bromide. 
G. Gloning and sequencing 

The identified differentially expressed bands were 
extracted from the agarose gel by using the PureLink 
Quick Gel Extraction Kit (Invitrogen). Each of these 
extracted DNA fragments was then individually cloned 



with the use of the TOPO TA Gloning Kit for 
Sequencing (Invitrogen) . Subsequently, the plasmid con- 
taining the inserted DNA fragment was extracted from 
clones of interest via PureLink Quick Plasmid Miniprep 
Kit (Invitrogen). The isolated cloned plasmids were 
then sequenced with the ABI 3730x/ DNA Analyser 
(Applied Biosystems). Finally, all the sequences obtained 
were analysed and matched for similarities with refer- 
ence to the BLAST programme under the NGBI 
database. 



RT-qPCR analysis 

A. Reverse transcription 

The total RNA isolated from 27 paired samples was 
reverse transcribed to first-strand cDNA, with the follow- 
ing protocol: 3 \ig of total RNA was added with 2 ]iL of 
0.5]ig/]iL oligo(dT) 12-18 (Invitrogen) and RNase-free 
water to a final volume of 9.5 ]iL. The reaction mixture 
was then incubated at 80 °G for 3 min, followed by chil- 
ling on ice for another 2 min. Next, 4 ]iL of 5X first- 
strand buffer (Invitrogen), 5 ]iL of 2 mM dNTP 
(Fermentas), 0.5 ]jiL of 40U/]jiL RNaseOUT recombin- 
ant RNase inhibitor (Invitrogen) and 1 \ih of 200 U/]jiL 
M-MLV reverse transcriptase (Invitrogen) were added to 
the mixture. Finally, the reaction mixture was incubated 
at 42 °G for 90 min, heated at 94° G for another 2 min 
and chilled on ice for 2 min. The synthesised first-strand 
cDNAwas stored under — 20°G until further usage. 

B. AAGx analysis 

The relative expression of identified DEGs in all paired 
colorectal tumours and control samples was determined 
via AAGx method. The RT-qPGR was performed in a sin- 
gleplex reaction containing 50 ng first-strand cDNA 
under universal thermal cycling conditions with the ABI 
7500 Fast Real-Time PGR System (Applied Biosystems). 
Both ACTB (Assay ID: Hs99999903_ml) and GAPDH 
(Assay ID: Hs99999905_ml) were used as reference 
genes and are commercially available as TaqMan 
Pre-designed Assays (Applied Biosystems). Prior to the 
analysis of gene expression, the amplification efficiency 
for all target and reference genes assays was measured 
by using the standard curve method with 2-log measure- 
ments. The amplification efficiency value of 90-110% 
was acceptable (Applied Biosystems). In this relative 
quantification method, the 2~^^^^ values obtained 
represented the fold change in gene expression of the 
colorectal tumours, which was normalised with both ref- 
erence genes, in relative to the calibrator (control 
sample) 

G. Statistical analysis 

The difference in the expression level between primary 
colorectal tumour and paired non-cancerous tissues was 
analysed by using Real-Time StatMiner software 
(Integromics). The distribution of the AGx values 
obtained for each DEGs within each GRG and control 
group were tested for normality via the Shapiro-Wilk 
test. Subsequently, the paired t test was performed to 
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assess the statistical significance of the observed differen- 
tial expression patterns. 

RESULTS 

DEGs between colorectal tumours and non-cancerous 
colonic tissues 

This preliminary study was conducted on paired samples 
pooled from four patients with CRC stages I-III. In 
ACP-based GeneFishing PGR, 20 sets of arbitrary AGP 
primers were used to randomly amplify gene products in 
both colorectal tumours and normal colonic samples. 
Upon visualisation on agarose gels, a total of 13 differen- 
tially expressed bands were observed by means of com- 
paring bands intensity between the tumouric and 
non-cancerous samples, as shown in figure 1. These 
bands were further sequenced for gene identification, 
and 16 DEGs were successfully reported. Among them, 
13 DEGs were overexpressed in colorectal tumours, 
while 3 DEGs were underexpressed, as listed in table 2. 



with stages I and II were grouped as early stage GRG, 
while the advanced stage GRG group comprised 14 
patients with stages III and IV. 

The analysis of RT-qPGR results was performed via 
Real-Time StatMiner software by importing the raw Gt 
data. The within-group correlation of these AGx values 
was then determined by calculating the median absolute 
deviation for all the samples within the same experimen- 
tal group. The biological samples that do not correlate 
well with other samples in the same group, were 
detected as group outliers and excluded from subse- 
quent analysis. Both ACTB and GAPDH were used for 
normalisation in computing the AGx (figure 2) and 
2-aact yalues by using the following formulas (table 4): 

^T(Target gene) ^(^(^Rgfgj-gnf-g gene) AGx 
^Cx(Sample) ~ ^Cx( Calibrator) = AAGx 

Relative fold change in expression (RQ) = 2~^^^^ 



Differential ability of the identified DEGs on early and 
advanced colorectal neoplasia 

Following the identification of DEGs, the gene 
sequences obtained were then used to design primers 
and TaqMan probes for RT-qPGR analysis by Applied 
Biosystems, as listed in table 3. In an attempt to assess 
the differential ability of identified DEGs on early and 
advanced colorectal adenocarcinomas, the recruited 
paired samples were further stratified into two groups 
according to the cancer stage. Among them, 13 patients 



The relative fold change in the mRNA expression 
level between the colorectal tumours and adjacent 
normal colonic mucosa was shown as the 2~^^^^ values. 
The statistical significance of the observed fold change 
in expression was determined by paired t test for all the 
DEGs. A p value of less than 0.05 is considered as statis- 
tically significant (table 4). 

In both early and advanced stage GRG groups, the 
expression of 4 of 16 DEGs was reported to be signifi- 
cantly differed between tumouric and non-cancerous 
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Figure 1 Differential banding patterns on 3% agarose gel postannealing control primer (ACP)-based PGR amplification between 
normal colon and colorectal tumour samples (N, normal sample; C, colorectal cancer sample). 
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Table 2 Sequence similarities and identification of DEGs ^ii^ 


Differentially 






Sequence 










expressed 






homology 


Accession 


UnlGene 




band 


DEG 


Identity 


(%) 


number 


number 


Description 


Overexpressed 
















A4.1 




Homo sstpions proteasorne (prosorne, 
macropain) 26S subunit, ATPase, 5 
(rblvlUo), mHNA 


(99%) 


NM_ 


_002805.4 


Ulo 7QQQ7 

MS./ yoo/ 


Involves in the ATP-dependent degradation of 
ubiquitinated proteins 


■ 


DEG2 


H. sapiens ubiquinol-cytochrome c 
reductase hinge protein (UQCRH), mRNA 


514/521 

/QQO/ \ 


NM_ 


_006004.2 


Hs.481571 


A component of the ubiquinol-cytochrome c 
reductase complex (complex III or cytochrome 
b-cl complex, which is part of the 
mitochondrial respiratory chain 


A4.2 


DEGS 


H. sapiens ribosonnal protein S23 
(RPS23), mRNA 


551/551 
(100%) 


NM_ 


_001 025.4 


Hs.527193 


A component of the 40S subunit of human 
ribosomes 


AD.1 


Ubbi4 


H. sapiens ribosomal protein L10 (RPL10), 
transcript variant 1, mRNA 


CCA /CC"7 

(99%) 


NM_ 


_006013.3 


LJo CO/1 Ar\A 


A component of the 60S subunit of human 
ribosomes 


AQ O 


UtLab 


H. sapiens actin related protein 2/3 
complex, subunit 2, 34 kDa (ARPC2), 
transcript variant 2, mRNA 


(100%) 


NM_ 


_005731 .2 


MS.iD^yoUo 


Involves in the regulation of actin 
polymerisation as an actin-binding component 
of the Arp2/3 complex, and mediates the 
formation of branched actin networks together 
with an activating nucleation-promoting factor 


r 


UtLa / 


H. sapiensT\W\P metallopeptidase inhibitor 
1 (TIMP1), mRNA 


(98%) 


NM_ 


_003254.2 


|_|o COO^iQO 


Irreversibly inactivates the metalloproteinase 
by binding to their catalytic zinc cofactor 


A10.1 




H. sapiens ATP synthase, H+ transporting, 
mitochondrial F1 complex, beta 
polypeptide (ATP5B), nuclear gene 
encoding mitochondrial protein, mRNA 


yi7/yiy 
(99%) 


NM_ 


_001 686.3 


rlS.40DOlO 


A subunit of mitochondrial ATP synthase that 
catalyses the synthesis of ATP by utilising an 
electrochemical gradient of protons across the 
inner membrane during oxidative 
phosphorylation 




DEG11 


H. sapiens chromosome 1 1 open reading 
Trame iu (uiiomuj, mHNA 


273/273 
(1 UU /o) 


NM_ 


_01 4206.3 


Hs.437779 


Unknown 


A13.3 
W 


DEG12 


H. sapiens mitochondrial ribosomal protein 
L24 (MRPL24), nuclear gene encoding 
mitochondrial protein, transcript variant 2, 
mRNA 


408/41 1 
(99%) 


NM_ 


_024540.3 


Hs.418233 


Involves in protein synthesis within the 
mitochondrion 


k A13.4 
1 


DEul3 


H. sapiens similar to UK/bW-OL.lD 

/I r\r^-i nnnoQ A -i Q\ 
(LUU1 UU^oo41 o) 


635/644 

/noo/ \ 
(yoTo) 


XM_ 


,002342023.1 




Unknown 


A18.1 


DEG14 


H. sap/ens family with sequence similarity 
96, member B (FAM96B), transcript 
variant 2, transcribed RNA 


486/487 
(99%) 


NR_ 


024525.1 


Hs.9825 


Involves in chromosome segregation as part 
of the mitotic spindle-associated MMXD 
complex 


|A20.1 


DEG15 


H. sapiens ribosomal protein L35 (RPL35), 
mRNA 


440/446 
(99%) 


NM_ 


_007209.3 


Hs. 182825 


A component of the 60S subunit of human 
ribosomes 


A20.2 


DEG16 






NM_ 


_001 01 2507.2 


Hs.486401 




Continued 



Open Access 



o 

(0 
0) 

o 



11 



8 E 



0) > 
O G) 
C O 

o- E ^ 
CO ^ ^ 



CM 

ll 



o 



o 

LU 
O 



(0 "O 
0) 

C (0 
0) (0 

£ Q. C 

S X 5 




03 rr ^ 
c/) O O 



c 

CO 

c 

CD 

O < 
CD I 

o ^ 

C/) CO 
O 

2 ^ 

^ CD 
<^ O 
CO ^ 

S=s CO 

ll 



CO 

1^ 

CO 

X 



Gi 

o 
o 
o 



CM 

o 



CO 

X 



CM 

00 
C^ 
lO 

o 
o 



CX) Co CD Co 
CM o^ 1- o^ 

^8^8 



T- CD ^ 



CM 
CO 
1^ 
O 

I 

cc 

X 



CD 




tissues. Remarkably, the combination of this panel of 
four genes is different among two groups. The RPL33, 
RPS23 and TIMPl genes were found to be overexpressed 
in early and advanced colorectal neoplasms (p<0.05) 
( figures 3 and 4) . It is interesting to note that the under- 
expression of ARPC2 gene (p<0.05) was only observed 
in early stage colorectal tumours (figure 3). However, 
the C6orfl73 gene was found to be overexpressed 
(p<0.05) in advanced colorectal adenocarcinomas, but 
not in early stage colorectal tumours ( figure 4) . 



DISCUSSION 

Our current study has revealed two distinctive four-gene 
signatures for both early and advanced stage colorectal 
adenocarcinomas. The early stage sporadic CRC was 
characterised by the overexpression of RPL35, RPS23 
and TIMPl genes, as well as underexpression of ARPC2 
gene. However, the advanced primary colorectal 
tumours were reported with overexpression of C6orfl73, 
RPL33, RPS23 and TIMPl genes. Although the relative 
fold change for ARPC2, RPL35 and RPS23 genes is 
below 2, the individual result does not affect the analysis 
as gene expression patterns of all four genes in combin- 
ation were proposed to distinguish between the early 
and advanced stage colorectal neoplasms. The potential 
involvement of these DEGs and their altered expression 
levels in CRC were further supported by previous 
researches. 

In fact, several proto-oncogenes and tumour suppres- 
sors are previously reported to regulate the ribosome 
production, that is, the RB}^ TP53}^ PTEN genes,^^ as 
well as the MYC gene family. It is suggested that the 
alterations in ribosome biogenesis might affect the trans- 
lation of genes that are involved in neoplastic transform- 
ation. In addition, the additional extra-ribosomal 
functions of the ribosomal proteins (r-proteins) in cellu- 
lar apoptosis, cellular proliferation, cellular transform- 
ation, genes transcription, mRNA translation, DNA 
repair and inflammation, might also trigger and support 
the neoplastic development.^^ Hence, the overexpres- 
sion of r-proteins-encoding genes observed in colorectal 
adenocarcinomas is not unexpected. Our current 
study has revealed the significant overexpression of two 
r-proteins that were not previously described in colorec- 
tal tumours, that is, the RPL35 and RPS23. The observed 
fold changes for the RPL35 and RPS23 mRNA levels 
were comparable between the early and advanced stage 
colorectal tumours in our sample cohort. This was in 
agreement with previous reports by Barnard et al^^ and 
Frigerio et al}^ where the changes in the mRNA expres- 
sion levels of the r-proteins were irrespective of the 
cancer stage. The hypothesis that the same ribosomal 
protein may contribute in different stages of cancer pro- 
gression with their hitherto unknown extra-ribosomal 
roles might provide an explanation to these 
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Table 3 Primers and TaqMan probes for relative quantification with Comparative Cj method. 





Primers sequence 


iac|ivian prooe sec|uence 
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FoHA/ard: 


5'-GGGCGTGTGCACAGAAG-3' 


5'-CTCGCAGGGCATACAT-3' 




Reverse: 


5'-AAGTCCTCCTGAGTGACATGGA-3' 




DEG2 


FoHA/ard: 


5'-GATGCTTACCGAATCCGGAGATC-3' 


5'-CCTCTTCCTCTTCCTCCTCC-3' 




Reverse: 


5'-GCATTGCTCTCTCACTGTTGTTAG-3' 




DEG3 


FoHA/ard: 


5'-CAACCGTCATTGGGTACAAAGG-3' 


5^-ATGGCAAGAAAATCAC-3^ 




Reverse: 


5'-TGTAAGGGTCCAGCTGATCAAGA-3^ 




DEG4 


FoHA/ard: 


5'-CGGCCAGGAAACTTGAACTTG-3^ 


5'-CAGGGCCTCAATCACA-3' 




Reverse: 


5'-CCGAGCTGCAGAACAAGGA-3' 




DEG5 


FoHA/ard: 


5'-CTGGTCGAATGAGGCACCTAAAA-3' 


5'-CATGGCTGAATCTGC-3' 




Reverse: 


5'-TGGGTTTAGGTGTTGTTCCTTCAC-3^ 




DEG6 


FoHA/ard: 


5'-AGATTAGCGGGATGAAAACGTCTT-3' 


5'-CCCCGTGATTG 1 1 1 IC-3' 




Reverse: 


5'-CGCCCAGATGCCGAGAAAA-3' 




DEG7 


FoHA/ard: 


5'-GGTAGTGATGTGCAAGAGTCCAT-3^ 


5'-CATTGCTGGAAAACTG-3' 




Reverse: 


5'-CCGCAGCGAGGAGTTTCT-3' 




DEG8 


FoHA/ard: 


5'-GAAGGAGACCATCAAAGGATTCCA-3' 


5^-ATTCACCTGCCAAAATC-3' 




Reverse: 


5'-GAAGGCCTGTTCTGGGAGATG-3' 




DEG9 


FoHA/ard: 


5'-GGCAGGGTGGTCCTGAGA-3' 


5'-CCTCTCTCCGCCCCGGACA-3' 




Reverse: 


5'-CCGCCATTGGCCTTAACTG-3' 




DEG11 


FoHA/ard: 


5'-CAGGTTTCAGTGAAGCCATCTG-3^ 


5'-CACCCAAGGGTAACAAC-3' 




Reverse: 


5'-GGGTTGGCATCTACGTGTGA-3' 




DEG12 


FoHA/ard: 


5'-CCAGGTCAAACTTGTGGATCCT-3' 


5'-ATGGACAGGAAACCCAC-3' 




Reverse: 


5'-GCTTCAGTAAATCTCCACTCGATCT-3^ 




DEG14 


FoHA/ard: 


5'-CCCGCTCCTTATCTGCAAGTT-3' 


5'-CATGCAGTGAACAAGC-3' 




Reverse: 


5'-TCAAGATGGACGTGCACATTACTC-3' 




DEG15 


FoHA/ard: 


5'-CGGCCTCCAAGCTCTCT-3' 


5'-CCGGACGACTCGGATCT-3' 




Reverse: 


5'-TGAGAACACGGGCAATGGATTT-3' 




DEG16 


FoHA/ard: 


5'-GGACTCTTCTGCTAATCGATGAACA-3^ 


5'-CAGATGGACCAATAAGTCA-3' 




Reverse: 


5'-GCCTCAACTTCGTCTGGAGAAAA-3' 




DEG, differentially expressed gene. 



2.95^ 
O 2.68- 
2.43- 
2.23- 
2.00- 
1.61- 



II 



Sample 
■ CRC 

DNormal 



ARPC2 \ CllorflO I FAM96B I PSMC5 I RPL35 I RPS23 I TIMP1 I 

ATP5B C6orf173 MRPL24 RPL10 RPL37 SLC25A1 UQCRH 
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ATP5B C6orf173 MRPL24 RPL10 



RPS23 y TIMP1 

SLC25A1 L 



Figure 2 Box plots showing ACj values of all colorectal tumours and normal colonic tissues in each early (A) and advanced (B) 
stage colorectal cancer (CRC) group. 
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Table 4 ACj mean, AACj, 2 and p values for all the differentially expressed genes (DEGs) in early and advanced stage 
colorectal cancer (CRC) groups 





Early stage CRC 








Advanced stage CRC 








i 


ACj mean 


ACj mean 








ACj mean 


ACj mean 








DEG 


(CRC) 


(normal) 


AACt 


2 -AACt 


p Value 


(CRC) 


(normal) 


AACt 


2 -AACt 


p Value 


AHPU2 


2.6854 


2.0664 


0.6190 


0.651 1 


0.0282 


2.7240 


2.3300 


0.3940 


0.7610 


0.2424 


AT Pod 


1 .5846 


1 .2702 


0.3144 


0.8042 


0.3524 


1 .9558 


1 .3838 


0.5720 


O.oJZJ 


0.1484 


UiiOniO 


3.2897 


3.3639 


—0.0742 


1 .0528 


0.8333 


3.3281 


3.6709 


-0.3428 


1 .2682 


0.3710 


C6orf1 73 


6.1083 


7.1943 


-1 .0860 


2.1228 


0.0905 


5.9949 


7.9087 


-1 .9138 


3.7680 


0.0013 


FAM96d 


3.5602 


3.8955 


-0.3353 


1 .2616 


0.2935 


3.5276 


3.9920 


-0.4644 


1 .3797 


0.21 13 


MHPL24 


4.9171 


5.0839 


-0.1668 


1 .1226 


0.3564 


4.9728 


5.1467 


-0.1739 


1 .1281 


0.7001 


rOIVIL/O 


3.8232 


3.9617 


— u. 1 ooo 


1 1 nop 
1 . 1 uuo 


U.Do 1 c. 


3.7705 


3.8455 


— u.u/ ou 


1 .UOoH- 


u.ouH-o 


RPL10 


-0.7462 


-0.4853 


-0.2609 


1.1982 


0.4001 


-1.1576 


-0.5196 


-0.6380 


1 .5562 


0.0950 


RPL35 


-0.1926 


0.6222 


-0.8148 


1.7591 


0.0024 


0.1748 


0.8769 


-0.7021 


1.6269 


0.0372* 


RPL37 


-0.0059 


-0.1539 


0.1480 


0.9025 


0.8645 


0.2184 


0.7143 


-0.4959 


1.4102 


0.1537 


RPS23 


0.2176 


0.7739 


-0.5563 


1 .4705 


0.0310* 


0.0676 


0.9431 


-0.8755 


1 .8346 


0.0250* 


SLC25A1 


3.7514 


3.5430 


0.2084 


0.8655 


0.5721 


3.5565 


3.4428 


0.1137 


0.9242 


0.7991 


TIMP1 


2.9096 


4.3059 


-1.3963 


2.6323 


0.0440* 


2.3330 


3.8547 


-1.5217 


2.8713 


0.0062* 


UQCRH 


2.0087 


2.2216 


-0.2129 


1.1590 


0.4108 


2.3375 


2.4459 


-0.1084 


1.0780 


0.7808 



*p<0.05=statistically significant. 



Figure 3 Differential expression 
patterns of all the identified 
differentially expressed genes in 
early stage colorectal cancer 
group. 
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Figure 4 Differential expression 
patterns of all the identified 
differentially expressed genes in 
advanced stage colorectal cancer 
group. 
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However, our current study also demonstrated an over- 
expression of the TIMPl gene in both early and 
advanced stage primary colorectal tumours. This finding 
is supported by Zeng et al^^ where the overexpression of 
TIMPl was reported in all stages of primary colorectal 
tumours. Under normal physiological conditions, the 
proteolytic activities of MMPs are kept at bay by their 
natural inhibitors, the TIMPs.^^ Previous studies have 
reported the overexpression of MMPs in early and 
advanced stage colorectal tumours, as well as other 
cancer types,^^"^^ which is in accordance to their bio- 
logical roles. Hence, a similar scenario is expected for 
TIMPs and indeed, their suppressive role in tumour 
invasion and metastasis has been demonstrated in 
various cancer models. However, more recent studies 
have revealed a direct correlation between TIMPl 
expression and tumour aggressiveness in cancer, includ- 
ing CRC.^^ These findings, which are contradictory to 
its protease-inhibiting function, have suggested a pos- 
sible tumour-promoting role of TIMPl in tumorigenesis. 
It is postulated that the TIMPl exhibits the abilities to 
inhibit tumour cell apoptosis and promote tumour 



angiogenesis, as well as other growth-factor-like effects. 
In our current study, the observed comparable overex- 
pression of TIMPl in early and advanced stage sporadic 
colorectal neoplasms was in line with its MMP inhibitory 
and MMP-independent tumour-promoting activities. 

In cancer biology, the expression of mRNAs and pro- 
teins of the ARP2/3 complex is often studied due to its 
role in cell migration, which contributes to cancer inva- 
sion and metastasis if aberrantly regulated.^^ We have 
detected a significant underexpression of ARPC2 in our 
cohort of early stage primary colorectal tumours. 
Surprisingly, this finding is contradictory with the role 
played by ARPC2 in cancer invasion and metastasis the- 
oretically. Previously, Kaneda et al has reported the 
decreased expression of all the seven genes encoding 
the subunits of ARP2/3 complex in human gastric 
cancers. Among them, the Arp2, ARPC2 and ARPC3 
showed the most prominent reduction in their expres- 
sion levels.^^ The exact mechanism underlying this 
observation still remains unknown, but the epigenetic 
alteration might potentially provide an explanation for 
it. For instance, promoter hypermethylation that causes 
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gene silencing is responsible for the reduced expression 
of ARPCl in human gastric cancer. Similarly, the epi- 
genetic study might also offer a clue for the underex- 
pression of ARPC2 in colorectal neoplasms. 

CSorfl 73, which is also known as CUG2 or CENP-W, is a 
novel oncogene that has been found to be upregulated 
in many human cancer tissues. Its high expression level 
is profoundly reported in tumours of the ovary liver, 
lung, pancreas, breast, colon, rectum and stomach. The 
CENP-W is a new member of the constitutive 
centromere-associated network, which specifically inter- 
acts with the CENP-T and plays an important role in 
mitosis. In our current study, the CENP-W is over- 
expressed in advanced colorectal adenocarcinomas. This 
finding correlates to its function in kinetochore assem- 
bly, where its aberrant expression might lead to abnor- 
mal cell division and aneuploidy in cancer. In our 
study, the overexpression of CENP-W was observed in 
early and advanced cohort of colorectal neoplasms but 
only statistically significant in the latter group. Given the 
fact that aneuploidy is constantly associated with a 
greater proportion of advanced CRC cases, the aberrant 
expression of CENP-W might potentially relate to a 
poorer prognosis of CRC.^^ 

In conclusion, we have characterised two distinctive 
gene expression patterns, which comprise the ARPC2, 
C6orfl73, RPL33, RPS23 and TIMPl genes, for the strati- 
fication of primary colorectal adenocarcinomas among 
Malaysian CRC patients. It was postulated that the actin 
cytoskeleton might play an important role in determining 
the dysplastic cell morphology during the early develop- 
ment of CRC, while the aberration in the assembly of 
functional kinetochore might be crucial for the aneu- 
ploidy of the advanced stage colorectal tumours. 
Nevertheless, the findings of this study were considered 
preliminary owing to the relatively small sample size. The 
main reason for this is the lack of a designated Tissue 
Bank in our institution. Moreover, the lack of patient with 
CRC volunteers and our stringent criteria for patient 
selection have also limited the availability of suitable spe- 
cimens within the short sample collection period. 

However, our identified mRNA expression patterns 
specific for early and advanced stage colorectal tumours 
are still convincing with our stringent sample selection 
criteria, high specificity primers and probes as well as 
reliable statistical analysis. In future, the validation of 
these DEGs should be performed on a larger set of clin- 
ical samples, and extensive inter-laboratory testing of 
their differential abilities on each CRC stage is also 
desired. In addition, we should also integrate other 
imaging and histological information to complement 
our identified gene expression patterns, which then 
hold promises for better stratification of tumours. 
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